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Abstract 

We study linearized perturbations of Myers-Perry black holes in d = 7, with two of the three 
angular momenta set to be equal, and show that instabilities always appear before extremality. 
Analogous results are expected for all higher odd d. We determine numerically the stationary 
perturbations that mark the onset of instability for the modes that preserve the isometries of 
the background. The onset is continuously connected between the previously studied sectors 
of solutions with a single angular momentum and solutions with all angular momenta equal. 
This shows that the near-extremality instabilities are of the same nature as the ultraspinning 
instability of d > 6 singly-spinning solutions, for which the angular momentum is unbounded. 
Our results raise the question of whether there are any extremal Myers-Perry black holes which 
are stable in d > 6. 
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1 Introduction 

In recent years, we have acquired a greater understanding of the phase space of higher-dimensional 
black holes. The main motivation for this research is that string theory predicts the existence 
of higher spacetime dimensions. Other scenarios with extra dimensions have been proposed in a 
phenomenological approach, independently of string theory. It is therefore crucial to understand 
what kinds of objects exist in higher- dimensional spacetimes. Black holes are the simplest gravita- 
tional objects, and they are essential to the study of fundamental theories. They are also the main 
candidates to signal the possible existence of extra dimensions in high-energy collisions. 

If higher- dimensional black holes can be produced in our universe, which are the relevant solu- 
tions? In four spacetime dimensions, the answer is simple: there is a unique vacuum solution, the 
Kerr black hole, which is (expected to be) stable. However, uniqueness fails in higher dimensions 
pp. The higher-dimensional version of the Kerr black hole, the Myers- Perry (MP) solution, has 
long been known [2] . It was the discovery of the black ring by Emparan and Reall [3] which showed 
that several rotating solutions may exist with the same asymptotic charges (mass and angular mo- 
menta), and that different topologies of the event horizon are possible. Another important feature 
is the possibility of black hole solutions with disconnected horizons, such as the black saturn [3] and 
concentric rings [5j|6], among others; this is believed to be impossible in four dimensions (in vac- 
uum). Actually, only the MP solution is known exactly for d > 5. The new solutions were obtained 
in d = 5, where powerful analytical techniques are available, but they are expected to exist in any 
higher dimensions. They have in fact been constructed approximately in a certain regime of high 
rotation [7], in what has recently been put into a systematic formalism - the blackfold approach 



An obvious question is whether these solutions are stable. There is one property of higher- 
dimensional black holes which is fundamental to this question: the event horizon can have very 
different length scales. In four dimensions, the horizon of the Kerr black hole is a deformed sphere, 
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and the deformation increases with the rotation. However, that deformation can never be too 
pronounced, because the rotation is bounded by extremality. This type of bound does not exist 
for some higher-dimensional solutions, so that regimes exist where there is a hierarchy between 
the scales set by the angular momenta and by the mass. When the mass scale is much larger, 
one expects a similar behaviour as in four dimensions. When an angular momentum scale is much 
larger, new features arise which are characteristic of higher dimensions. 

The approximation of large rotation, which is the basis of the blackfold approach mentioned 
above, was explored by Emparan and Myers to argue that singly-spinning MP black holes are 
unstable if rotating too rapidly [UJ. Singly-spinning MP solutions have the isometry group R x 
U(l) x SO(d — 3), which is an enhancement from the isometry group R x U(l) n of generic MP black 
holes (n = \_(d — l)/2j is the number of independently defined angular momenta). The argument 
for instability is that, for very high rotation, the solution acquires a disk- like horizon and resembles 
a black brane locally. It has been shown by Gregory and Laflamme [12] that black branes are 
unstable, so these ultraspinning MP black holes should become unstable too. Ref. [TT] also pointed 
out that, for perturbations which do not break the U(l) symmetry (which is the angle associated 
to the rotation of the black hole), the onset of the instability should be a stationary mode which 
signals the bifurcation to a new black hole family. This idea was elaborated on Ref. [7], which 
conjectured that there is a harmonic- type structure of ultraspinning instabilities, and that the new 
black hole solutions appearing at the associated bifurcations interpolate between the MP family 
and the black ring, the black saturn, concentric rings, etc.. 

Refs. |13|ll4j provided strong evidence for this picture of the phase space of singly-spinning black 
holes. The critical rotations and the stationary onset modes were determined numerically. More- 
over, those modes were found to be consistent with the expected horizon shape of the new families, 
if they are to connect to the black ring, black saturn, etc.. Ref. [13] put forward a conjecture about 
the appearance of these instabilities beyond the singly-spinning sector. It says that instabilities 
whose onset is a stationary mode can only appear if the black hole has two local thermodynamic 
instabilities. Note that this is a necessary, but not sufficient condition for the appearance of such an 
instability. We refer to this conjecture as the ultraspinning conjecture. It seems strange to relate 
classical and local thermodynamic stability for non-extended objects (as opposed to the Gubser- 
Mitra conjecture for black branes [IS])- The point is that these instabilities are part of the same 
harmonic structure. The two lowest harmonics are associated with the asymptotic charges and 
with the local thermodynamic instabilities, while the higher harmonics are associated with modes 
which cannot change the asymptotic charges, and are associated with the classical instabilities. 

The goal of the ultraspinning conjecture is to provide a guidance for the appearance of instabil- 
ities. The great advantage is that it does not rely on hierarchies between the mass and the angular 
momenta scales, telling us something about the intermediate regime where the scales are of the 
same order. The conjecture was the motivation for the work of Ref. [16], which studied the stability 
of cohomogeneity-1 MP black holes. This sector consists of odd-dimensional solutions which have 
all n angular momenta turned on, but equal. There is an enhancement of symmetry to the isometry 
group R x U(N + 1), where d = 2N + 3; this leads to a line element which depends non-trivially 
on the radial coordinate only, i.e. which has cohomogeneity-1 [17] . The regime of large rotation, 
in the sense that the angular momentum is much larger than the scale set by the mass, does not 
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exist in this sector, because the rotation is bounded by extremality. However, the ultraspinning 
conjecture allows for a small parameter region close to extremality where an instability is possible, 
in d > 7. The d = 5 solution is stable [18] . 

An instability in the cohomogeneity-1 sector was indeed found in d = 9 [16J. Moreover, due 
to the symmetry of these solutions, it was possible to determine the growth rate of the unstable 
modes, which increases with rotation, while the works on the singly-spinning sector could only find 
the stationary onset. This was an explicit confirmation that we are dealing with actual instabilities. 
Ref. [16J failed to find that the instability appears already in d = 7. It is so close to extremality, 
(<^inst — Jext)/Jext ~ 10~ 5 , that it did not show up in the original analysis. It was revealed in a 
reanalysis performed by one of us (JES), by request from the authors of Ref. [19], where it is argued 
that the extremal black hole is unstable. The instability extends to a small parameter region close 
to extremality, as the numerics confirmed. Ref. |19j makes use of an extension of the Newman- 
Penrose formalism to higher dimensions [20], and studies the stability of extremal solutions by 
analysing their near-horizon geometry. It is conjectured there that the existence of an instability 
of the near-horizon geometry implies the existence of an instability of the full black hole geometry, 
if the near-horizon perturbations satisfy a certain symmetry. 

The story seems satisfactory so far, and the ultraspinning conjecture has been shown to hold 
also in the singly-spinning asymptotically AdS case |21j . where a large rotation limit analogous to 
the asymptotically flat case exists |22[ I23|. However, an important question remains. Why should 
the instability in the cohomogeneity-1 sector be related to the instability in the singly-spinning 
sector? While in the latter case there is the convincing argument of Emparan and Myers, no 
such clear geometric understanding of the instability exists for the cohomogeneity-1 sector. One 
can see the ultraspinning conjecture as extending the notion of brane-like behaviour, and thus the 
argument of Emparan and Myers. However, it is important to see explicitly how these two sectors 
fit together in the stability picture of the MP family. This provides the motivation for the present 
paper, where we will study a MP sub-family which includes the cohomogeneity-1 sector and the 
singly-spinning sector as particular cases. 

We will consider MP black holes in an odd number of spacetime dimensions, which have all but 
one of the n angular momenta equal, so that we have effectively two unequal spins. The isometry 
group of this MP sub-family is 1 x U{1) x U(N) (recall that d = 2N + 3). We particularize our 
study to d = 7, and restrict ourselves, for simplicity, to perturbations which preserve the isometries 
of the background solutions. As a consequence, the perturbed geometries will also have two unequal 
spins. This restriction still allows us to find the missing link between the singly-spinning and the 
cohomogeneity-1 sectors. One of the main results of this paper is that the onset of the ultraspinning 
instability is continuously connected between those two sectors in the MP parameter space. This 
shows that the instabilities are indeed of the same nature. 

The most important outcome of our analysis, however, goes beyond the initial motivation. We 
find that all near-extremal solutions studied here are unstable. This raises the question of whether 
there is any extremal MP black hole which is stable in d > 6, i.e. when the parameter space admits 
unbounded sectors. 

Before proceeding, let us point out that fastly rotating MP black holes can also be unstable 
against bar-mode perturbations that break the rotation-generating symmetry of the background. 
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This is the case at least for singly-spinning MP black holes \24\ I25j . We emphasize that the 
ultraspinning conjecture only applies to instabilities whose onset is a stationary mode, which is not 
the case for the bar-mode perturbations. Therefore, the bar-mode instability avoids the conjectured 
bound, and even occurs in d = 5, but it is not associated with a bifurcation to a new family of 
stationary black holes. It would be important to determine how this type of instability manifests 
itself beyong the singly-spinning sector. In d = 5, Ref. |18j found no evidence of instability for 
equal spins. 

This paper is organised as follows. In Section[2j we present the MP sub-family with two unequal 
spins, and show how the ultraspinning conjecture applies to it. In Section [3j we describe the linear 
perturbation problem that is solved numerically. We conclude in section |4| with the discussion of 
the results. 



2 Myers-Perry black holes with two unequal spins 
2.1 Solution 

The Myers-Perry (MP) black hole [2] is an asymptotically flat solution to the vacuum Einstein 
equations, which generalizes the Kerr solution to higher dimensions. This family of black holes can 
be parameterized by the asymptotic charges of the spacetime, the mass M and the n = [(d — l)/2j 
angular momenta Jf, or, equivalently, by the mass-radius tm and the rotation parameters aj, to be 
defined below. 

A generic MP solution has isometry group R x U(l) n , where R corresponds to time translations 
and each U(l) describes one of the n rotational isometries. However, there is an enhancement of 
symmetry when some of the angular momenta coincide. We will be interested in the sub-family 
of MP black holes that has n — 1 equal rotation parameters, denoted by a, while the remaining 
rotation parameter, denoted by b, is independent. In odd-dimensional spacetimes, which we will 
focus on, such solutions have the isometry group R x £7(1) x U(N), where d = 2N + 3. This family 
interpolates between the singly-spinning MP black hole and the equal angular momenta MP black 
hole. The former is obtained by setting a = 0, and has isometry group R x U(l) x SO(d — 3). 
The latter is obtained by setting a = b, and (for odd d) has isometry group R x U(N + 1); due 
to this enhancement of symmetry, this solution is cohomogeneity-1, i.e. it depends non-trivially on 
the radial coordinate only. 

For concreteness, we study here the MP sub-family with two unequal spins in d = 7. (In d = 5, 
there is no instability of the type which we will consider.) We then have a\ = a2 = a and 03 = b. 
The line element of this solution, which can be derived from [26], is 
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where we defined the functions 



x ( ,)-^±^!±g>--Ju, rw— -° 2) , (2.2) 



and 



d4 Pl = | (<i(9 2 + sin 2 (9 ciy? 2 ) , Ai = - cos dp, (2.3) 

are, respectively, the line element of CPi = S 2 and the associated Kahler 1-form Ai (related to the 
Kahler 2-form Ji by dA\ = 2Ji). The time t and radial r coordinates have the standard range, 
— oo < t < +oo, < r < oo, and the angular coordinates v, ipi, ip2, 0, 4> range in the intervals 

min{|a|, |6|} < v < max{|a|, |6|}, < tpx , V>2 < 2vr , < 6 < it , < 99 < 2vr. (2.4) 

The mass M and the two unequal angular momenta J^ 1 ) and j(^ 2 ) are given by 

Note that a d = 7 spacetime has three asymptotically defined angular momenta (n = 3). However, 
since two of them are equal, it is more convenient to define from the line element (2.1) two angular 
momenta, with respect to d^ )1 and d^ 2 . 

The event horizon is located at r = r+, where r + denotes the largest real root of X(r) = 0, and 
it is a Killing horizon of £ = dt + difa + 0^f 2 ^<9^ 2 , where the angular velocities of the horizon 
are given by 

)<<'i) _ 0" -Q- Q W> 2 ) 
a(r2 +a 2)' 

Finally, the Hawking temperature Tjj and the Bekenstein-Hawking entropy S of the solution are 
given by 

r_+ ( r\-g 2 | 1 \ _ ^r A M r+ 

i ^"2vrl v r 2 .(r 2 +a2) + r2+62 y )' ^ " 4 ' ^ 

In terms of these quantities, the first law of thermodynamics takes the form 

dM = T H ds + n^dj^ + n { ^ 2) dj^ . (2.8) 

Since we defined the angular velocities/momenta with respect to and d^ 2 , this form of the first 
law does not allow for variations of parameters in the MP family which go beyond the two unequal 
spins sub- family; that is, this expression restricts the variations to be within the line element (2.1). 
This is important for us, because we will make use of the thermodynamics to indicate the possibility 
of classical instabilities of the black hole, and we will consider only perturbations preserving the 
equality between two of the three angular momenta. 

We presented here the general two unequal spins solution. Since our goal is to connect previous 
work on the singly-spinning sector and on the equal spins sector, we will analyze those particular 
limits in detail. The coordinate transformations 

R 2 1 b 

t = T-ail>, r=^/p 2 -a 2 , v = a + (b - a) — — ^ , ^i = x*, ^2 = -^, (2.9) 

1 + H A 2 a 
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followed by the limit b — > a take the MP solution (2.1) into the equal angular momenta MP black 
hole ITT]: 



do 2 

ds 2 = -f(p) g(p) dT 2 + jjrjr + P 2 h(p) [dip + A 2 - Q(p) dT] 2 + p 2 ds 2 p2 , where 
di? 2 i? 2 i? 2 i? 2 

ds " 2 = (TTW + (T+W + Al)2 + 1+^ ds ^ pl ' A2 = (d * + Al) ' 

Here, ds 2 p2 is the Fubini-Study metric on CP 2 , while A 2 is the associated Kahler 1-form. In this 
solution, surfaces of constant T and p have the geometry of a homogeneously squashed S 5 , written 
as an S 1 fibred over CP 2 . The fibre is parameterized by the coordinate tp, which has period 2ir. 
The largest real root of /(p) = gives the location of the event horizon, p = p + , which is a Killing 
horizon of £ = dt + ^Ind^ with £Ih = aj p\. The angular momentum of this black hole is bounded, 
being maximal for a regular extremal solution (the temperature vanishes but not the entropy). 
On the other hand, the coordinate transformations 

t = T-b<j>, r = r, v = bcos6, -d} 1 = ^ + -6, ip 2 = 4> , (2.11) 

a 



followed by the limit a — > 0, take the MP solution (2.1) into the singly-spinning MP black hole [2]: 



ds 2 = - ^ (dT - b sin 2 9 #) 2 + ^ [(r 2 + b 2 )dcf> - b dT] 2 + | dr 2 + £ ^ 2 + r 2 cos 2 dfl| 3 , 

2-j Zj _A 



A(r) = r 2 + 6 2 -^f , £(r, 0) = r 2 + 6 2 cos 2 , d0 2 3 = (dtf + A x ) 2 + ds 2 pl . (2.12) 

Note that dO| 3 is the line element of a unit <S 3 , here written as a Hopf fibration of S 1 over CP 1 = S 2 . 
The largest real root of A(r) = gives the location of the event horizon, r = r+, which is a Killing 
horizon of £ = + Vtud^ with $7// = ^rrp-- For the singly-spinning MP black hole (in d > 5), 
there is no bound on the angular momentum. 

Let us also mention the case 6 = 0. There is an extremality bound on the rotation parameter 
a of such a solution, but the limit is singular (both the temperature and the entropy vanish). In 
summary, the rotation parameters of the two unequal spins MP sub-family are bounded by regular 
extremal solutions, except for a = (where no extremal limit exists) and for 6 = (where an 
extremal limit exists, but is singular). 
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2.2 Thermodynamic zero-modes and the ultraspinning regime 

The ultraspinning conjecture proposed in Refs. |X3|, [16] says that classical instabilities of vacuum 
black holes whose onset is a stationary and axisymmetric mode can only appear in a parameter space 
region (called the ultraspinning region) determined by the existence of two local thermodynamic 
instabilities. By axisymmetric we mean that the symmetry generated by the Killing vector field 
VLiWii is preserved; here, the m^s denote the rotational Killing vector fields and the fVs denote 
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the associated angular velocities of the horizon. If axisymmetry is broken, then stationarity is not 
possible, by virtue of the rigidity theorem (27J EE]- I* 1 this section, we review this conjecture, which 
is essential to interpret our results. 

The condition for local thermodynamic stability is the positivity of the thermodynamic Hessian 

-S°fl = -P^, x Q = (M,J l ), (2.13) 

where S is the black hole entropy, M its mass and Jj denote all the possible angular momenta. 
Using only the Smarr relation and the first law, it was shown in Ref. [16] that —S a p possesses at 
least one negative eigenvalue for any asymptotically flat vacuum black hole. Hence all such black 
holes are locally thermodynamically unstable. The ultraspinning region is defined as the parameter 
space region where there are at least two negative eigenvalues of the thermodynamic Hessian. 
Let us elucidate the conjecture by considering MP black holes. For small angular momenta, 



there is a single negative eigenvalue of the Hessian (2.13), which is the one continuously connected to 



the negative specific heat of the Schwarzschild black hole. However, for d > 5, additional eigenvalues 
vanish and then become negative as the rotation increases; the vanishing of the first of these signals 
the ultraspinning surface. So the ultraspinning surface defines a "small angular momenta" region 
in the MP parameter space, where a single eigenvalue is negative, and no instability (whose onset is 
a stationary mode) is allowed. Beyond the ultraspinning surface, we have the ultraspinning region, 
where instabilities are allowed, but not required. 

Whenever the thermodynamic Hessian has a vanishing eigenvalue, the corresponding eigenvector 
(SM, 5Ji) represents a perturbation of the black hole (within the MP family, in our example) which 
changes its mass and angular momenta, but preserves its temperature and angular velocities |16j . 
This type of perturbation is associated with the £ = and £ = 1 "harmonics" of the black hole; 
these contribute to the asymptotic charges. What then corresponds to the £ > 2 harmonics? If 
such zero-mode perturbations were allowed, they would change the solution but not the asymptotic 
charges, as they decay too quickly with large radius. They would then indicate a bifurcation to a 
new black hole family which can have the same asymptotic charges of the unperturbed solution. In 
the MP case, such new families of black holes where conjectured to exist (in d > 5) and to provide 
a connection in parameter space between the MP family and other families of black holes, such as 
black rings, black saturns, etc. [7J[29]. Moreover, the MP black hole is expected to become unstable 
for rotations larger than that critical value of bifurcation [7J. The ultraspinning conjecture reflects 
the expectation that a zero- mode associated with the £ > 2 harmonics, and thus with the onset of 
the classical instability, appears only for rotations larger than the £>1 harmonics. 

In fact, there is a precise sense in which the harmonics £ > 2 can acquire a negative eigenvalue, 



just as the harmonics £ = and £ = 1. The spectrum of the thermodynamic Hessian (2.13) 
is generalized by the spectrum of the Euclidean action for (stationary and axisymmetric) off- 
shell perturbations of the black hole. This spectrum determines the one-loop correction to the 
gravitational partition function; see |30| 131] for the study of these corrections. That negative 
eigenvalues of the thermodynamic Hessian imply negative modes of the Euclidean action has been 
shown in Ref. [32], in the static case, and in Refs. [331 EE], i n the general case; the proof is based on a 



preceding construction [Ml [35]. These £ = and I = 1 type modes, predicted by the Hessian (2.13) 
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are called thermodynamic negative modes. Additional negative modes in the partition function are 
non-thermodynamic negative modes, and those are the ones associated to classical instabilities of 
the black hole. 

Ref. |32j pointed out that the eigenvalue equation for the Euclidean negative modes of a vacuum 
black hole, which we present in the next section, is the same as the equation for the threshold sta- 
tionary mode of a Gregory-Laflamme instability of a vacuum black brane, which trivially extends 
that black hole along extra dimensions. The negative eigenvalue of the Euclidean negative mode 
corresponds to (minus) the 'mass-squared" coming from the dimensional reduction along the brane 
directions. Therefore, a local thermodynamic instability of the black hole, i.e. a negative eigenvalue 
of the thermodynamic Hessian, implies the existence of a classical instability of the black brane. 
The simplest example is the original Gregory-Laflamme instability of a Schwarzschild black brane, 
which is related to the negative apecific heat of the Schwarzschild black hole. This is what under- 
lies the Gubser-Mitra conjecture [15J in the vacuum case. That conjecture says that black branes 
which have a translational invariance are classically unstable if and only if they are locally thermo- 
dynamically unstable. However, the preceding discussion shows that the Gubser-Mitra conjecture 
can be refined: there should be a distinct Gregory-Laflamme type instability for each negative 
mode (thermodynamic or not). Therefore, we will be finding not only thermodynamic and classical 
instabilities of the MP black hole, but also the associated classical instabilities of the MP black 
branes. 

The important point for our present work is that we do not expect an instability of the black 
hole (whose onset is a stationary mode) within the ultraspinning surface, in the MP parameter 
space. In the asymptocally flat case, the ultraspinning surface is defined as the locus of the first 
zero-mode of the thermodynamic Hessian — S a p, as we mentioned. However, such a zero-mode is 
also identified by a reduced thermodynamic Hessian, H^ = — (^ gfgj. 1 = — Sij, as shown in [16] . 
We can thus make use of the reduced Hessian analysis for the general MP black hole done in |14j . 

We will consider here only perturbations preserving the equality between two of the angular 
momenta. The MP ultraspinning region is defined having in mind perturbations which can change 
the angular momenta in an arbitrary manner, i.e. considering all the eigenvalues of Hij. Therefore, 
we can find a more stringent ultraspinning surface (i.e. a smaller ultraspinning region) for our 
perturbations if we use the angular momenta defined with respect to and d^ 2 , so that Hij is 
2x2 matrix, rather than 3x3. Specializing it to the case d = 7, and a\ = 02 = a and 03 = b that 
interests us, we find that the parameter space locus of the ultraspinning surface is: 



unequal spins sub-family is presented (ignore the coloured dots for now). Inside this surface (i.e. 
to the left of the black curve in Figure [TJ, — S a p has only one negative mode. There is a second 
eingenvalue of — S a p that is positive inside the ultraspinning surface, but then vanishes and becomes 
negative as the rotation increases and the ultraspinning surface is crossed (i.e. when we cross the 
black curve in Figure [T] from left to right). 





This ultraspinning surface is the black curve in Figure [TJ where the parameter space of the two 
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The ultraspinning region is bounded by the ultraspinning surface and by the extremality curve 
where the black hole temperature vanishes. The parameter space locus of the extremality curve is 



M = v 2! /f + 1 ' ( 2 - 15 ) 

r M V b z 



and is described by the blue curve in Figure [TJ 



3 Ultraspinning perturbations of MP black holes 
3.1 Perturbations and eigenvalue problem 

We will consider perturbation modes h^ v that preserve the M x £7(1) x U (2) isometries of the d = 7 
MP solutions with two unequal spins (2.1 ). This means that we will look for stationary modes which 
preserve the equality between two of the three angular momenta of the background geometry. The 
most general perturbed line element with these properties is described by the ansatz: 
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(3.1) 



where {6a, 6f3, 5x, 6/11,6^2, 6^3, 6fJH, 6001,6002,6003} are small quantities that describe our perturba- 
tions, and they are only functions of (r, v). 

We choose to work in the traceless-transverse (TT) gauge, 



= , and Vhp, 



0, 



for which the linearized Einstein equations are 



9 n p <r u 



o. 



(3.2) 



(3.3) 



where is the Lichnerowicz operator. Actually, following |X3|, I16j. we will solve the more general 
eigenvalue problem 

{A L h)^ = -k 2 h^ . (3.4) 



This problem arises when one considers the stability of the MP black string, ds 2 tr - mg = g^ u dx^dx" + 



dz 2 , under Gregory-Laflamme-type perturbations, ds 2 



string 



ds 2 tri + e lkcZ h /lu (x)dx^dx lJ , where 



is the metric of the MP black hole ( |2.1[ ) Q The same problem also arises when one considers the 



1 The label k c , rather than k, means that this is the critical wavenumber of a Gregory-Laflamme perturbation. 
The perturbation decays with time for jfcj > k c , and grows with time for |fe| < k c . The stationary onset is |fe| = k c . 
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quadratic quantum corrections to the gravitational partition function in the saddle point approxi- 
mation |31| (see |36j for the application to the Kerr-AdS black hole, where the numerical method 
used here was first applied to black hole perturbations). 



Equation (3.4) describes a coupled system of partial differential equations (PDEs) that we can 



solve only numerically. We look for solutions of (3.4) instead of (3.3) because Mathematica has 



powerful built-in routines to solve generalized eigenvalue problems. The strategy to find zero- modes 
is then to look for solutions of (3.4), and vary the rotation parameters (a, b) until we obtain k c = 0. 



An added value of this strategy is that we will be looking not only for instabilities of the MP black 
hole but, in addition, to instabilities of the associated black strings. Indeed, solutions with k c ^ 
signal new kinds of Gregory-Laflamme instabilities and inhomogeneous phases of ultraspinning 
black strings. On the other hand, when k c = 0, we will be dealing either with the thermodynamic 
zero-mode (whose location in parameter space - the ultraspinning surface - we predicted in the 
last section), or with the onsets of the ultraspinning instabilities of the MP black hole. 



Notice that the ansatz in (3.1) is the most general one that respects the isometries of the 



background MP black hole (2.1) and is preserved under diffeormorphisms that depend only on 
(r, v). Ultimately, this is necessary and sufficient to guarantee that (3.4) leads to a closed system 
of equations. 



3.2 Boundary conditions 

In this section, we will discuss in detail the boundary conditions that we need to impose on the 



metric perturbations in order to solve (3.4). In the present situation, we have to specify boundary 

\a\ and v = \b\. 
To find the associated 



conditions at the horizon, r = r+, at asymptotic infinity, r — > oo, and at v 
The metric perturbations must be regular on the future event horizon 7i~ 
boundary conditions we first introduce the ingoing Eddington-Finkelstein coordinates (V, ifti, ife): 



dt = dV 



r 2 + a 2 



X(r) 



dr . 



dtpi = difii + 



aX(r) 



dr . 



dip 2 = dip2 



b(r 2 +a 2 ) 
r 2 X(r) 



dr . 



(3.5) 



In these coordinates, the background geometry (2.1) is manifestly regular on % . Regularity of 
the perturbed geometry is guaranteed if the components hav of the perturbation in the Eddington- 
Finkelstein coordinates are regular. Translated to the perturbation notation introduced in the 



ansatz (3.1), the boundary conditions at the horizon are then 



5a, 5/3, 5ni, 5fi 2 , , ^3, #M4> Scjx, 5bj 2 = 0(1) 



Su 3 , d~x = 0(r - r+) . 



(3.6) 



Quite importantly, note that perturbations obeying these boundary conditions preserve the tem- 
perature and the angular velocities of the background black hole. This is the case because we 
have imposed the regularity of the perturbation h liV dx iL dx v separately, i.e. seen as a symmetric 
2-tensor on the unperturbed background, rather than imposing only the regularity of the perturbed 
geometry. 

At spatial infinity, r — > oo, we are interested in boundary conditions that preserve the asymp- 
totic flatness of the spacetime. In this asymptotic region, the Lichnerowicz eigenvalue equations 
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(3.3) reduce simply to Dh^ ~ k 2 h/j, u . Its regular solutions decay as 



Ku\ e~ kcr ^0, (3.7) 

for some constant a > that depends on the particular metric component. Therefore, for k c ^ 0, 
asymptotic flatness is preserved since the perturbations vanish exponentially. (Note that, in our 
numerical method, we find modes that approach as close as desired k c = 0, but without ever 
reaching it exactly.) 

The strategy to discuss the boundary conditions at v = \a\ and v = \b\ is similar to the one 
used for the event horizon. We first zoom the background geometry near these critical points in 
coordinates where the geometry is manifestly regular. The boundary conditions for the metric 
perturbations can then be determined by demanding that the components hu, v are regular in the 
above coordinates where the background metric is explicitly regular. We leave the details of this 
analysis to Appendix^ and quote the final conclusion here. 

The boundary conditions at v = \a\ are 

5a, 5(3, 5 hi, 5fi2, 5uj\, 5oj2, 5u)3, 5\ = 0(1) j <fy*3 — 5fi4, 5^a — 5fii = 0(v — a) . (3.8) 

On the other hand, the boundary conditions at v = \b\ are 

5a, 5/3, 5[ii, 5n2, 5ui2, 5\ = 0(1) , (3.9) 

a 2 

5fi3 ~ Sm = 0(v - b) , Sui - -= = 5uo 2 = 0(v - b) 5u 3 = Ofv - b) . 

r J, _|_ a z 



As a verification of our regularity analysis, we have checked that the boundary conditions (3.6) 



(3.10) are consistent both with the Lichnerowicz eigenvalue equations (3.3) and with the TT gauge 



conditions (3.2). Indeed, the first term in the series expansion of the eigenvalue Lichnerowicz 



equations vanishes after we impose the aforementioned boundary conditions. Moreover, the first 



term of a series expansion of the TT gauge conditions is also consistent with (3.6)-(3.10). 



3.3 Strategy to find the unstable modes 

The strategy to solve the Lichnerowicz eigenvalue problem (|3.4|) for the ten metric perturbations 



described in (3.1), namely {6a, 5/3, Sx, 8fJ,i, SfJ,2, 8fJ>3, 5fj,±, 5u\, 5uj%, 5us}, subject to the TT gauge 



conditions (3.2) is as follows. 



We first introduce the orthogonal tetrad basis {e^}: 



e (°)=dt+^ — ^3-(^i+Ai), e«=dr, e^=dv, e ( 3 ) = dt + ° \ + °^ {d^ x + A x ) , 
a z — b z ar — b z 

( 1 ( r 2 + g?) C„2 _ v 2\ 1 \ 

e( 5 ) = rvab -^dt- V , 2 J L 2 J (#1 + A x ) + -or #2 h e« = #*> & dx* , (3.10) 



a (a 2 — b 2 ) r 2 v 2 a?b 



where is the vielbein of the CP 1 = S 2 manifold. In this particular tetrad basis, the Lichnerowicz 
eingenvalue equations simplify considerably. 
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The map between the metric perturbations h^ v in the coordinate basis and the perturbations h a \, 
in the tetrad basis follows straightforwardly from h a b = e, ^e, h ^hn V . The ten metric perturbations 



{5a, 5(3, 5x, 5fxi,5/j,2, 5^l 3 , <5//4, 5u\, 5lo 2 , 5oj 3 } described in the ansatz (3.1 ) then map to the ten tetrad 
components {h 00 , h 02 , h 03 , tin, ti\ 2 , h 22 , h 33 , h 3i , /i 44 , h 55 }. 

The TT gauge conditions (3.2) allow us to eliminate {/ioch ^03> ^34} in terms of the other seven 
functions {Hq 2 , hn, h\ 2 , h 22 , h 33 , /144, /155} and their first derivatives. Making these substitutions in 



the full set of the perturbation equations (3.4), we find that only seven equations remain of second 
order. Explicitly, these equations are 



-klh\ 2 , 



(A L h) 



22 



-k 2 h 22 , 



(A L h) 02 = -k c h 02 , (A L h) u = -k c h u , (A L h) l2 
(A L h) 33 = -k 2 c h 33 , (A L h) u = -k 2 h M , (A L h) 55 = -k 2 c h 55 , (3.11) 

and they constitute our final set of equations to be solved. A non-trivial consistency check of our 



procedure is to verify that the final equations (3.11) imply that the three remaining equations 



which are of third order, are automatically satisfied. 



We use spectral methods to solve numerically the final eigenvalue problem (3.11). Less com- 



putational power is required in the implementation of the spectral method if all functions obey 
Dirichlet boundary conditions on all boundaries. Thus, we introduce new independent functions 
{qi} related to the independent metric functions by 



qi(r,v) 
q 3 (r,v) 

qi{r,v) 
q 5 (r,v) 



a)(v - b) (l - r -±) 2 h u (r,v) , q 2 (r,v) = ±-(v - a)(v - 6) (l - ^) h l2 (r,v) , 

a) 2 ( W -6) 2 (l-^) h 22 (r,v), 

J)2 h 33 (r,v) - Y{v) 2 h 22 {r,v) 



(r 2 + v 



( r 2 +v 2f Y {v) 



-{v - a)(v -b)(l- T -yj hu(r,v) , 

, _ b) 2 ( 1 _ ! + \ (a 2 ~ b 2 ) 2 h 55 {r, v) - a 2 Y(v) 2 h 22 (r, v) 
V ' V r J (a 2 -b 2 ) 2 Y(v) 



q 7 (r,v) = h 02 (r,v) . 

(3.12) 



The Dirichlet boundary conditions for the qiS are equivalent to the boundary conditions (3.6) 
(3.7), (3.8) and ( 3. lOh PI It is also convenient to work with coordinates that range in a unit interval 



We thus introduce the new coordinates 

a — v 



y 



1 



r 



x 



with a = min{|a|, |6|} and f3 = max{|a|, , 



(3.13) 



a — (3 

so that < y < 1 and < x < 1. 

2 The reader might ask the reason for the particular combinations in the definitions of q± and qe, in contrast with 
the simpler definitions of the other qi's. These combinations maximize the information on the boundary conditions. 



To illustrate this take (3.10 1. One has 6/13 — 6(11 = 0(v — b). So our fundamental variables should be {6fii, 6/13 — 6/11}, 
instead of considering only the leading behaviour of {6/11, 6/13}; otherwise the information that 6/13 — 6/ii — is lost. 
A similar reasoning applies to the variables involved in the last relation of (3.8 1. The combinations chosen for the 



definition of 54 and qe account for this issue, in the notation of the tetrad metric components. 
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12 3 4 

b/r + 



Figure 1: Phase diagram of MP solutions with two unequal spins in d = 7. The axes are the 
rotation parameters a and b, normalized to the horizon radius r+. The colour code is the number 
of negative modes found numerically: red 1, blue 2, violet 3, green 4, brown 5, yellow 6. The blue 
line represents extremality, the black line between red and blue represents the ultraspinning surface, 
and the line between the blue and violet regions is the onset of the ultraspinning instability. The 
scattered white dots are numerical faults. 



4 Results and discussion 



We determined numerically the negative modes in the Lichnerowicz eigenvalue problem (3.4) for the 
MP sub-family with two unequal spins in d = 7. Our perturbations /i™ are stationary and preserve 
the U(l) x U(2) spatial isometry of the background solution, which implies that they preserve the 
equality between two of the three MP spins. In Figure [TJ we present the number of negative modes 
for a grid of points on the parameter space. The axes represent the rotation parameters, a and 
b, normalized to the horizon radius r+. The blue line that binds the parameter space represents 
extremality. The colour coding is the following: red denotes the existence of 1 negative mode, blue 
2, violet 3, green 4, brown 5, and yellow 6. 

For small rotation parameters, in the red region, there is a single negative mode. This is the 
negative mode that any asymptotically flat vacuum black hole possesses [IB] : for a = b = 0, it 
is the negative mode of the Schwarzschild black hole [31], which is associated with its negative 
specific heat. In the blue region, an additional negative mode exists. This negative mode is also 
thermodynamic in nature. It appears along the black line, which represents the ultraspinning 



surface, discussed in Section 2.2 This surface, given by the expression (2.14), is defined by the 



appearance of a second negative eigenvalue of the thermodynamic Hessian (2.13). Therefore, the 



fact that it precisely separates the red and blue regions indicates an excelent agreement between 
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our expectations and the numerical results. 

The most interesting line in Figure [l] is the separation between the blue and the violet regions J^] 
Across that line, a third negative mode appears, one which is not related to the thermodynamic 
Hessian (2.13). This mode is the first stationary onset of the ultraspinning instability (within our 
symmetry restrictions of keeping two unequal spins) , and it marks also a bifurcation to a new family 
of black holes. The prediction of the ultraspinning conjecture is verified, as this line only appears 
for rotation larger than the ultraspinning surface. The lines between the regions violet/green, 
green/brown and brown/yellow (and others which would appear had we zoomed in the appropriate 
region, or extended our analysis to higher b/r + ) signal new instabilities/bifurcations. Let us also 
point out that each of those lines (and also the black one, the ultraspinning surface) corresponds 
to the appearance of a new Gregory-Laflamme-type instability of the MP black branes. 

One of the main results of this paper is that the blue/violet line, i.e. the onset of the first 
ultraspinning instability, continuously connects the singly-spinning case (a = 0) and the equal 
spins case (a = b). These sectors have been studied before in Refs. [13] and [Ml [19], respectively, 
and the numerics matches precisely with the previous results. The question which remained was 
whether the instabilities in the singly-spinning sector and in the equal spins sector were of the 
same nature. We recall that there is a geometric argument that explains the instability in the 
singly-spinning case |llj . namely that for high rotation the black hole resembles locally a black 
brane, which suffers from the Gregory-Laflamme instability. However, no analogous argument has 
been proposed for the instability in the equal spinning sector, for which the rotation is bounded 
by extremality. We have shown here that the instability in the equal spins sector is indeed of the 
same nature as the originally proposed instability in the singly-spinning sector. Their onsets are 
continuously connected in the MP parameter spacej^] 

What is possibly our main result is that the blue/ violet onset line seems to occur everywhere 
before extremality. This would mean that we can always find unstable black holes in a small 
enough neighbourhood of extremality, and that the extremal black holes are themselves unstable. 
This raises several questions. 

We can follow numerically only a finite number of onset lines, e.g. blue/violet, violet/green, 
green/brown, etc.. There are then three possibilities which are difficult to check numerically. The 
first is that all the onset lines squeeze in between the ultraspinning surface and extremality, beeing 
closer to extremality for higher number of associated negative modes. This would mean that the 
extremal black holes would be severely unstable solutions, since they would suffer from an infinite 
number of distinct linear instabilities. In the second and third possibilities, there is an infinite 
number of onset lines which terminate at some point along the extremality line, but in a different 



3 The reason for choosing the horizon radius as a normalization is that the "distance" between the blue region 
and extremality is larger, say for a = b. Had we normalized the grid with respect to the mass-radius tm, it would 
be impossible to see the region of interest for equal spins without zooming in substantially. On the other hand, this 
normalization makes the parameter space seem unbounded for b = 0, because r+ — > in the singular extremal limit. 

4 Note that Refs. [161 119] consider general axisymmetric perturbations in the equal spins sector, without our 
symmetry restriction of keeping the spatial isometry (7(1) x (7(2). The onset of the ultraspinning instability found 
here is a single perturbation, while in those works there was a large degeneracy of the harmonics of CP 2 from which 
the perturbations were constructed. This means that several "onset surfaces" will meet at the equal spins point in 
the full MP parameter space, and we have studied only one of them, which was sufficient for our purposes. 
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manner. (In fact, we can't be sure that this does not occur for the blue/violet line too, for small 
b/r + .) So the second possibility is that the onsets can occur also on the extremal line, i.e. there is 
a regular zero- mode for some extremal solutions, which would connect them to new extremal black 
holes. This is not possible for perturbations of a scalar field on an extremal black hole background, 
since the scalar field cannot be regular on the near-horizon geometry [37] ■ The third possibility 
is that the same is true for gravitational perturbations, in which case the onset lines terminate 
at extremality but the zero-mode is not regular for the extremal solution, so that the bifurcation 
does not extend to extremality. For instance, in the scalar field condensation analysis of Ref. |38j, 
the zero-mode line touches the extremality line tangentially, but at that point one does not have a 
hairy extremal black hole. 

The most relevant question raised by the instability of extremal solutions, in the MP sub-family 
studied here, is whether this extends to MP extremal solutions beyond this sub-family. Recall that, 
in d = 5, no instability was found on the equal spins sector [18]. However, only in d > 6 are there 
regimes with unbounded angular momenta, so the parameter space is qualitatively different. So 
the question is: are there any extremal Myers-Perry black holes which are stable in d > 6? 

This is the first work which studies the stability of MP solutions with two independent angular 
momenta turned on. As we emphasized, we focus on axisymmetric perturbations. It would be 
interesting to know how perturbations breaking that symmetry, such as the ones studied in Refs. |24} 
I25| . constrain the stable parameter space. Moreover, we still lack a good understanding of the 
stable parameter regions of other higher-dimensional solutions, for instance black rings. The main 
problems in higher-dimensional vacuum black holes remain open: (i) Is a classification of higher- 
dimensional black holes possible? (ii) Which MP solutions are stable? (iii) Are there any stable 
solutions apart from "not-so-rapidly rotating" MP black holes? 
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A Boundary conditions at v = \a\ and v = \b\ 



In this Appendix, we give the details of the regularity analysis at v = a and v = b (take a, b > 



for simplicity), which leads to the boundary conditions (3.8) and (3.10). 

Let us start with the analysis of v = a. For v « a, we can write Y(v) = Y'(a)(v — a) + 0[(v — a) 2 ] 
with Y'(a) = 2{b 2 — a?) /a. To zoom the geometry in this region, we introduce the new variable p: 
v = a + p 2 / 4. We introduce also the manifestly regular 1-forms E p = pdp and E^ 1 = p 2 (dtpi +Ai). 
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The geometry near v = a then reads 



ds z 



r 2 + a 2 
X(r) 



r 2 a ( r 

dr 2 + 



2(b 2 - a 2 ) 



Js 2 

as ball 4 



2(6 2 



air 2 + a 2 



E^dt + p 2 



(b 2 



,2\2 



r 2 _|_ Q 2 



dt + 



2(b 2 - a 2 ) 



E^ 



+ r 2 a% 2 



#2 

a 2 b 2r 2 a 2 {b 2 - a 2 ) 



a 2 {r 2 + a 2 ) 2 
EH 



dr 



(A.l) 



r 2 + a 2 



dt 



r I _J_ a z 



where we used ds 2 aUi = dp 2 + p 2 [(dtpi + Ai) 2 + ds 2 p i]. This line element is manifestly regular 
as p — > 0. The boundary conditions for the metric perturbations can now be determined by 
demanding that h^dx^dx" is a regular symmetric 2-tensor, i.e. that it has regular components 
when expressed in the above coordinates, where the background metric is regular. Near p = 
(v = a), the perturbation reads 



h^y dx^ dx v 



r 2 + a 2 



X(r 
f2ab 2 (r 2 + a 2 ) 



{2813 dr 2 - 5 X E p dr) + 

2 „2 



a(r 2 +a 2 ) 



b 2 -a 



2E^dt + * \, p 2 dt 2 
air* + a ) 



2 

2„4l2 



+ 2r z a% z 5p 



(dpi dp + 5p, 3 p (dtpi + Ai) + 5//4 ds cp i) 

o 



#2 . 

21 + 



r 2 + a 2 
2r 2 a 2 (6 2 - a 2 ; 



+ 



Sun E^ 1 - 4a 2 b 2 5u 2 dt 



b 2 



dvp2 . r 2 + a 2 



+ 



2 b 2r 2 a 2 (b 2 -a 2 ] 



E^ 



dt 



E^ 1 

2 



r 2 a 2 



r 2 a 2 



2aX(r) E?h 
r z + a z o 

2X(r) 



r 2 + a 2 



8a [dt + 



2(6 2 - a 2 ) 
E^ 



E^ 1 ] <Jw 3 + 2di ( E^ 1 



p dt \ 0W3 



a(r 2 + a 2 ) 



2(6 2 - a 2 



(A.2) 

Regularity of the perturbed geometry as v — > a then requires the boundary conditions (3.8). Note 



that the last condition in (3.8) guarantees regularity of the last contribution in between curved 



brackets in the first line of (A.2); see the contribution ds 2 aUi in the background geometry ( A.l[ ). 

Finally, consider v = b. In the neighborhood of v = b, we can write Y(v) = Y'(b)(v — b) + 0[(v — 
b) 2 ] with Y'(b) = 2(a 2 — b 2 )/b. The new variable p, defined through v = b — p 2 /A, allows us to zoom 
the geometry in this region. Introduce the new coordinates dt = dt + bdij)2 and dt/i\ = dip% — - dt/j2- 
Consider also the regular 1-forms E p = pdp and E^ 2 = p 2 dip2- The background geometry, in the 
vicinity of v = b, is then 



ds z 



r 2 + b 2 , 2 b(r 2 + b 2 ) 
dr 1 + 



X(r 
+2 



ar 2 + a 2 



b r 2 + b 2 



+r z [(dip! +Ai) 



2(6 2 - a 2 ) 
(#1 + Ai 
a 



ds 



q2 (y,2 _|_ q2\2 
balls + 52 ( r 2 +fe 2)2 ^ ( (#1 + A l) 



a(r 2 + a z 



-dt 



a(r 2 + a 2 ) 



dt) E^ 2 



dt 



X(r) 
r 2 + b 2 



dt + a(dipi + Ai) 



+ (r + a )ds C¥l 
2 



(A.3) 



where we used ds 2 all = dp 2 + fPdtp^- This is a manifestly regular geometry as p — > 0. Requiring 
that the components hp V are regular when expressed in the these coordinates fixes the boundary 
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conditions at v = b (p = 0). In the neighborhood of v = b, the perturbation reads 



h dr^ dr u « - ; 

n^ux ux 2{b 2 - a 2 ) 



5ji\ dp 2 + 5p, 3 p 2 I dip2 + 



ar 2 + a 2 
b r 2 + b 2 



+ Ai 



b 2 -a 2 / 2 
dt 



b(r 2 + b 2 ) 



4 r 

5bJ; 



(r 2 + a 2 )(di)i + Ai) -adt 



4aX(r) 
r 2 + b 2 



-dip 2 + #1 + A] 



(r 2 + a 2 )5wi ( -dip 2 + rf^i + Ai ] - a5u 2 (bdip 2 - dt) 



+ 25p 2 



dt + a(dipi + Ai 



(r 2 + a 2 )(dVi + Ai) - adt 



2(r 2 + a 2 ) ~ r ~ ~ 

^ 6uj 3 u9 ._ 7 ,_ 9 . , (&#2 - dt) (r 2 + a 2 )(dVi + Ai) - adt 



+ 2 



6 2 r 2 (r 2 + 6 2 ) 
r 2 + b 2 
X(r) 



(<5/3 dr 2 + 5x drdv) + 2(r 2 + a 2 )(5// 4 ds^ - 5a 



2X(r) 
r 2 + b 2 



dt + a(d^i + Ai) 



(A.4) 

It follows that the perturbed geometry is regular as v — > b if and only if the boundary conditions 
(3.10) are obeyed. Note that in the second line of (3.10), the first, second and third conditions 
guarantee the regularity of the first, second and third lines of (A.4), respectively. 



References 

[1] R. Emparan and H. S. Reall, "Black Holes in Higher Dimensions," Living Rev. Rel. 11 (2008) 
6, |arXiv:0801 .34711 [hep-th]]. 

[2] R. C. Myers and M. J. Perry, "Black Holes In Higher Dimensional Space-Times," Annals Phys. 
172, 304 (1986). 

[3] R. Emparan and H. S. Reall, "A rotating black ring in five dimensions," Phys. Rev. Lett. 88 
(2002) 101101 |arXiv:hep-th/0110260|. 



[4] H. Elvang and P. Figueras, "Black Saturn," JHEP 0705 (2007) 050 |a rXiv:hep-th/0701035] . 

[5] H. Iguchi and T. Mishima, "Black di-ring and infinite nonuniqueness," Phys. Rev. D 75 (2007) 
064018 [Erratum-ibid. D 78 (2008) 069903] |arXiv:hep-th/0701043| . 

[6] J. Evslin and C. Krishnan, "The Black Di-Ring: An Inverse Scattering Construction," Class. 
Quant. Grav. 26 (2009) 125018 [arXiv: 0706 .12311 [hep-th]]. 

[7] R. Emparan, T. Harmark, V. Niarchos, N. A. Obers and M. J. Rodriguez, "The Phase 
Structure of Higher-Dimensional Black Rings and Black Holes," JHEP 0710 (2007) 110 
[arXiv:0708.2T8T1 [hep-th]]. 

[8] R. Emparan, T. Harmark, V. Niarchos and N. A. Obers, "Blackfolds," Phys. Rev. Lett. 102 
(2009) 191301 | ;arXiv:0902.0427| [hep-th]]. 

[9] R. Emparan, T. Harmark, V. Niarchos and N. A. Obers, "Essentials of Blackfold Dynamics," 
JHEP 1003 (2010) 063 [arXiv:0910.160T1 [hep-th]]. 



17 



[10] R. Emparan, T. Harmark, V. Niarchos and N. A. Obers, "New Horizons for Black Holes and 
Branes," JHEP 1004 (2010) 046 [arXiv:0912.2352l [hep-th]]. 

[11] R. Emparan and R. C. Myers, "Instability of ultra-spinning black holes," JHEP 0309 (2003) 
025 |arXiv:hep-th/0308056| . 

[12] R. Gregory and R. Laflamme, "Black strings and p-branes are unstable," Phys. Rev. Lett. 70, 
2837 (1993) [arXiv:hep-th/9301052] . 

[13] O. J. C. Dias, P. Figueras, R. Monteiro, J. E. Santos and R. Emparan, "Instability and 
new phases of higher-dimensional rotating black holes," Phys. Rev. D 80 (2009) 111701 
[arXiv:0907.2248l [hep-th]]. 

[14] O. J. C. Dias, P. Figueras, R. Monteiro and J. E. Santos, "Ultraspinning instability of rotating 



black holes," Phys. Rev. D 82 (2010) 104025 |arXiv: 1006. 1904 [hep-th]]. 



[15] S. S. Gubser and I. Mitra, "Instability of charged black holes in anti-de Sitter space," arXiv:hep- 
|th/0009126[ S. S. Gubser and I. Mitra, "The evolution of unstable black holes in anti-de Sitter 



space," JHEP 0108 (2001) 018 arXiv:hep-th/0011127l|. 



[16] O. J. C. Dias, P. Figueras, R. Monteiro, H. S. Reall and J. E. Santos, "An instability of 
higher-dimensional rotating black holes," JHEP 1005 (2010) 076 |arXiv:1001.4527l . 

[17] H. K. Kunduri, J. Lucietti and H. S. Reall, "Gravitational perturbations of higher dimensional 
rotating black holes: Tensor Perturbations," Phys. Rev. D 74 (2006) 084021. |arXiv:hep^| 
th/0606076] . 



[18] K. Murata, J. Soda, "Stability of Five-dimensional Myers-Perry Black Holes with Equal An- 
gular Momenta," Prog. Theor. Phys. 120 (2008) 561-579. [ arXiv:0803.1371| [hep-th]]. 

[19] M. Durkee, H. S. Reall, "Perturbations of near-horizon geometries and instabilities of Myers- 
Perry black holes," Phys. Rev. D83 (2011) 104044. [arXiv:1012.4805] [hep-th]]. 

[20] M. Durkee, H. S. Reall, "Perturbations of higher-dimensional spacetimes," Class. Quant. Grav. 
28 (2011) 035011. [arXiv:1009.0"015l [gr-qc]]. 

[21] O. J. C. Dias, P. Figueras, R. Monteiro, J. E. Santos, "Ultraspinning instability of anti-de 
Sitter black holes," JHEP 1012 (2010) 067. |arXiv:1011.0996l [hep-th]]. 



[22] M. M. Caldarelli, R. Emparan, M. J. Rodriguez, JHEP 0811 (2008) 011. |ar Xiv:0806.1954 
[hep-th]]. 



[23] J. Armas, N. A. Obers, Phys. Rev. D83 (2011) 084039. | arXiv:1012.5 081 [hep-th]]. 



[24] M. Shibata, H. Yoshino, "Nonaxisymmetric instability of rapidly rotating black hole in five 



dimensions," Phys. Rev. D81 (2010) 021501. arXiv:0912.3606 [gr-qc]]. 



18 



[25] M. Shibata, H. Yoshino, "Bar-mode instability of rapidly spinning black hole in higher di- 
mensions: Numerical simulation in general relativity," Phys. Rev. D81 (2010) 104035. 
[arXiv: 1004.4970| [gr-qc]]. 

[26] T. Oota, Y. Yasui, "Separability of Gravitational Perturbation in Generalized Kerr-NUT-de 
Sitter Spacetime," Int. J. Mod. Phys. A25 (2010) 3055-3094. |arXiv:0812.l"623| [hep-th]]. 

[27] S. Hollands, A. Ishibashi, R. M. Wald, Commun. Math. Phys. 271 (2007) 699-722. (gr^l 
|qc/0605T06j . 



V. Moncrief, J. Isenberg, Class. Quant. Grav. 25 (2008) 195015. (arXiv:0805.T45l] [gr-qc]]. 



[29] R. Emparan, P. Figueras, "Multi-black rings and the phase diagram of higher-dimensional 
black holes," JHEP 1011 (2010) 022. [arXiv: 1008.3243] [hep-th]]. 

[30] G. W. Gibbons, M. J. Perry, "Quantizing Gravitational Instantons," Nucl. Phys. B146 (1978) 
90. 

[31] D. J. Gross, M. J. Perry and L. G. Yaffe, "Instability Of Flat Space At Finite Temperature," 
Phys. Rev. D 25 (1982) 330. 

[32] H. S. Reall, "Classical and thermodynamic stability of black branes," Phys. Rev. D 64 (2001) 
044005 | jarXiv:hep-th/0104071] . 

[33] R. Monteiro, M. J. Perry and J. E. Santos, "Thermodynamic instability of rotating black 
holes," Phys. Rev. D 80 (2009) 024041 |arXiv:0903.3256l [gr-qc]]. 

[34] B. F. Whiting and J. W. . York, "Action Principle and Partition Function for the Gravitational 
Field in Black Hole Topologies," Phys. Rev. Lett. 61 (1988) 1336. 

[35] T. Prestidge, "Dynamic and thermodynamic stability and negative modes in Schwarzschild- 



anti-de Sitter," Phys. Rev. D 61 (2000) 084002 arXiv:hep-th/9907163|. 



[36] R. Monteiro, M. J. Perry and J. E. Santos, "Semiclassical instabilities of Kerr-AdS black 
holes," Phys. Rev. D 81 (2010) 024001 [arXiv:0905.2"334l [gr-qc]]. 



[37] J. Fernandez-Gracia, B. Fiol, JHEP 0911 (2009) 054. arXiv:0906.2353 [hep-th]]. 



[38] O. J. C. Dias, R. Monteiro, H. S. Reall, J. E. Santos, JHEP 1011 (2010) 036. [arXiv: 1007.3745 
[hep-th]]. 



19 



